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ABSTRACT 

The observed power spectrum of redshifted 21cm fluctuations is known to be sen- 
sitive to the astrophysical properties of the galaxies that drove reionization. Thus, 
detailed measurements of the 21cm power spectrum and its evolution could lead to 
measurements of the properties of early galaxies that are otherwise inaccessible. In this 
paper, we study the effect of mass and redshift dependent escape fractions of ionizing 
radiation on the ability of forthcoming experiments to constrain galaxy formation via 
the redshifted 21cm power spectrum. We use a model for reionization which com- 
bines the hierarchical galaxy formation model GALFORM implemented within the 
Millennium-II dark matter simulation, with a semi-numerical scheme to describe the 
resulting ionization structure. Using this model we show that the structure and distri- 
bution of ionised regions at fixed neutral fraction, and hence the slope and amplitude 
of the 21 cm power spectrum, is dependent on the variation of ionising photon escape 
fraction with galaxy mass and redshift. However, we find that the influence of the 
unknown escape fraction and its evolution is smaller than the dominant astrophysi- 
cal effect provided by SNe feedback strength in high redshift galaxies. The unknown 
escape fraction of ionizing radiation from galaxies is therefore unlikely to prevent mea- 
surement of the properties of high redshift star formation using observations of the 
21cm power spectrum. 
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1 INTRODUCTION 

In recent years a great deal of theoretical attention has fo- 
cused on modeling the effect of galaxies on the reionization 
of the IGM. In large modern simulations, the most common 
approach is to begin with an N-body simulation to generate 



a distribution of dark matter halos (e.g. Ciardi et al.||2003 



Sokasian et al.|2003|[riiev et al.|2007||Zahn et al.|2007||Trac 
fc Cen||2007| |Shin et al.||2008| flliev et al.||2008| [Trac et ~ 
2008[ ). A prescription is then used to relate dark matter halo 
mass to ionizing luminosity. Following this step, radiative 
transfer methods (most commonly ray-tracing algorithms) 
are employed to model the generation of ionized structure 
on large scales. The radiative transfer is normally run with 
lower resolution than the N-body code for computational 
efficiency. 



An important outcome from the large cosmological vol- 
umes attained by modern numerical simulations has been 
the prediction of 21 cm signals that will be observable us- 



ing forthcoming low frequency arrays (e.g. |Mellema et al 



2006 



Lidz et al.||2008[ ). |Lidz et al.| ( |2008[ ) argue that first 



generation low frequency radio telescopes like the Murchi- 
son Widefield ArrajJ^] and the Low Frequency Arrajj^] will 
have sufficient sensitivity to measure the redshift evolution 
in the slope and amplitude of the 21 cm power spectrum. 
One of the main limitations in modelling of reionization is 
the physics of the ionizing sources. Most studies have used 
very simple prescriptions to assign ionizing luminosities to 
dark matter halos. It has been shown that it is then possible 
to constrain the parameters for these simple prescriptions 
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(e.g. Barkana 2009). However, an important open question 
is the degree to which the important astrophysics governing 
formation and evolution of high redshift galaxies is accessi- 
ble via observations of the 21cm power spectrum. 

Several studies have previously addressed the issue of 
realistic modelling of high redshift galaxies and their role in 



reionization. For example, Raicevic et al. (20111 (see also 
Benson et al. 2006| Lacey et al. 20111 used the hierarchi- 



cal galaxy formation code, GALFORM (Cole et al. 2000 



Baugh et al. 2005; Bower ct al. 20061, implemented within 



Monte-Carlo merger trees to evaluate the ionizing photon 
budget. They found that although galaxies should produce 
sufficient ionizing photons to complete reionization, most 
of the galaxies responsible would be below the detection 
threshold of current surveys. However, these studies were 
restricted to global evolution of the ionized fraction, and 
did not address the ionization structure of the IGM. Most 



recently, |Kim et al.| ([2013} have combined the GALFORM 
model implemented within the Millennium-II dark matter 
simulation, with a semi-numerical scheme to describe the 
resulting ionization structure. In this paper we extend this 
model to consider an escape fraction for ionising photons 
that depends on redshift and host halo mass. 

The escape fraction of ionizing photons from their host 
galaxies is one of the most important unknowns for the 
reionization history. Observational estimates show a broad 
range of escape fraction values from a few percent in the local 



Unvierse (Hurwitz et al. 111997 Bland-Hawthorn & Maloney 



1999 Putman et al. 20031, to a few percent or a few tens 



of percent at redshift z ~ 1-3 (Inoue et al. 2006 Shapley 



et al. 2006 Chen et al. 2007 Siana et al. 2007). However, 



there are no observational constraints on the escape fraction 
during the epoch of reionization. Theoretically, the escape 
fraction is predicted to span a very broad range from 0.01 to 



1 (Wood & Loeb 2000; Ricotti & Shull||2000| |Benson et al. 


2002 


|Fujita et a .||2003| |Sokasian et al.| 2003| |Ciardi et al. 


2003 


Wyithe & Loeb|2007||Wise & Cen|2009||Raicevic et al. 


2011 


Kuhlen & Faucher-Giguere|2012 \. 



In this paper we incorporate a variable escape fraction 
for ionising photons within the model of Kim et ah] ( |2013[ ), 
and predict the redshifted 21cm power spectrum for the 
resulting reionization histories. We begin in § [2] and § [3] 
by describing the implementation of GALFORM, and our 
method for modelling the ionization structure. Then, in § [4] 
we present variant models which lead to different reioniza- 
tion histories, and present a discussion of the dependence of 
the 21cm power spectrum on the escape fraction models in 
§ [5] We finish with some conclusions in § [6] 



2 THE MODEL 

In this section we summarise the theoretical galaxy forma- 



tion modelling based on (Kim et al. 20131 that is used in 



our analysis in order to describe the new features for this 
paper. 

2.1 The GALFORM galaxy formation model 

We compute the formation and evolution of galaxy prop- 
erties using the hierarchical galaxy formation model GAL- 
FORM. GALFORM includes a range of processes that are 



thought to be important for galaxy formation from dark 
matter merger histories to baryonic physics such as cool- 
ing, star formation, feedback processes, and chemical en- 
richment. A comprehensive overview of GALFORM can be 



found in Cole et al.| (|2000[) , with an updated discussion in 
the review by Baugh (20061. Feedback processes (including 



AGN, SNe and photo-ionization) during galaxy formation 
play very important roles in shaping the luminosity func- 



tions predicted by GALFORM (Cole et al.||2000| 


Benson 


et al. 2002 


|Baugh et al.|2005||Bower et al.|2006||Kim et al. 


2011 12013 


|. In this paper, we use the version of GALFORM 



lation on which the model is based has a cosmology with 
fractional mass and dark energy densities of fi m = 0.25 , 
Sib = 0.045 and S7a=0.75, a dimensionless Hubble constant 
of h=0.73, and a power spectrum normalisation of 08=0.9. 
The particle mass in the simulation is 6.89xl0 6 h _1 M© and 
we detect halos down to 20 particles in the simulation box 
of side length L = 100/i _1 Mpc. 



SEMI-NUMERICAL SHEME TO 
CALCULATE THE EVOLUTION OF THE 
IONISED STRUCTURE 



Mesinger & Furlanetto (20071 introduced an approximate 



but efficient method for simulating the reionization process. 
In this paper we apply a semi-numerical technique to find 
the ionization structure resulting from GALFORM galaxies 
within the Millennium-II dark matter simulation, following 



the procedure described in Kim et al. (2013) 



We assume the number of photons produced by galaxies 
within a small volume (or cell) of the simulation, and which 
enter the IGM to participate in reionization to be 

iV 7 , C eU=/ esc / " ^Lyc,ceu(t) ttt, (1) 
JO 

where / csc is the escape fraction of ionizing photons pro- 
duced by stars in a galaxy. The total Lyman continuum lu- 
minosity of the iVceii galaxies within the cell expressed as the 
emission rate of ionizing photons (i.e. units of photons/s) 
is 



Ai 



Lyc,cc 



where 



N^ycAt) = 



hv 



(2) 



(3) 



L Ui i is the spectral energy distribution of galaxy i, and thresh 
is the Lyman-limit frequency (/ifthresh = 13.6 eV). Note that 
the number of photons produced per baryon in long-lived 
stars and stellar remnants depends on the IMF and metalic- 
ity (Z). We assume the total Lyman continuum luminosity in 
a cell at redshift Zj to be constant until the next snapshot at 
redshift Zi+i , and calculate the number of photons produced 
in the cell between Zi and Zj+i as iVLyc.ceiiCtzJ x (fs i+ i — 

We then calculate the ionization fraction within each 
cell according to 



N. 



7, cell 



(1 + FcWm.ceU 



(4) 
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where F c denotes the mean number of recombinations per 
hydrogen atom up to reionization and -/Vm.ceii is the number 
of neutral hydrogen atoms within a cell. The latter quantity 
is calculated as 



A^Hi.ceii = hm(Sr>M,cen + l)V cc ii, 



(•») 



where we assume that the overdensity of neutral hydrogen 
follows the dark matter (computed based on the Millennium- 
II simulation density field), rim is the mean comoving num- 
ber density of hydrogen atoms, and V ce n is the comoving 
volume of the cell. Self-reionization of a cell occurs when 
Qccii > 1. We divide the Millennium-II simulation box into 
256 3 cells, yielding cell side lengths of 0.3906/t _1 Mpc and 
comovmg volumes of 0.0596/i" 3 Mpc 3 . 

Theoretical prediction of the parameters F c and / eS c 
in Equations and Q from first principles is compli- 
cated, and their values are not known. The recombination 
parameter F c is related to the density of the IGM on small 
scales, while / esc depends on the details of the high redshift 
ISM. Previous work using GALFORM suggested the value 
(1 + F c )/fe SC ~ 10 ( |Benson et al.|2001||Raicevic et al.|20TT 
Kim et al.|2013 1 to fit observational constraints on reioniza- 



tion. Here, we vary the value of f esc , assuming F c =0.5, as 
a function of redshift or host dark matter halo mass to see 
the effect of different histories of reionization on the 21cm 
power spectrum. 

Based on equation ffl, individual cells can have Q C eii > 
1. On the other hand, cells with Q cc n < 1 may be ionized by 
photons produced in a neighbouring cell. In order to find the 
extent of ionized regions we therefore filter the Q cc ii field us- 
ing a sequence of real space top hat filters of radius R (with 
0.3906 < R < lOO/i^Mpc), producing one smoothed ion- 
ization field Qr per radius. At each point in the simulation 
box we find the largest R for which the filtered ionization 
field is greater than unity (i.e. ionized with Qr > 1). All 
points within the radius R around this point are considered 
ionized. This procedure forms the position dependent ion- 
ization fraction ^ Q ^ 1, which describes the ionization 
structure of the IGM during reionization. 

In this paper we restrict our attention to analyses that 
assume the spin temperature of hydrogen is coupled to the 
kinetic temperature of an IGM that has been heated well 
above the CMB temperature (i.e. T s 3> Tcmb). This condi- 
tion should hold during the later stages of the reionization 
era (z < 9, |Santos et aT1|2008| ). In this regime there is a 
proportionality between the ionization fraction and 21 cm 
intensity, and the 21 cm brightness temperature contrast 
may be written as 



AT(z) = T (z) [1-Q](1 + 5 E 



(6) 



where Tq(z) = 23.8 ( ]rf) 5 mK. We have ignored the con- 
tribution to the amplitude from velocity gradients, and as- 
sumed as before that the hydrogen overdensity follows the 
dark matter (<5dm,cc,ii). 



4 MODELLING THE HISTORIES OF 
REIONIZATION 

In this paper, we assume several different models for / eBC 
as defined in § [3] Specifically we model dependencies of the 
escape fraction with halo mass and redshift using 



Table 1. The values of selected parameters which are different 
in the models. The columns are as follows: (1) the name of the 
model, (2) Normarlization value to match < Xi >=0.5 at z=7.272, 
(3) index a describing redshift dependence, (4) index /3 describing 
halo mass dependence. 





A 


a 




Model 


0.5348 








Model I 


0.1488 


5 





Model II 


1.6791 


-5 





Model III 


0.3649 





5 


Model IV 


0.8966 





-5 



A x 



1 + z 



(logic (M ha io/h- 1 M )/lO) /3 J (7) 



where A is normalisation parameter, a describes the redshift 
dependence, and j3 describes the halo mass dependence. We 
force to f eB c~0 when the value of Eq.[7]is negative and / CS c = l 
when the value greater than unity. We also consider a fidu- 
cial model with / ea c=constant. We assume F c to equal 0.5 
(Barkana & Loeb||2001| iFurlanetto & Oh||2005h. To facili- 



tate comparison, we assume that A is normalized so that 
the ionization fraction at z=7.272 is (a;i)=0.5 for all models. 
In Table. [I] we show the parameters describing the 5 models 
used in this paper. These models are discussed below. 



4.1 Redshift dependence modelling 

There is uncertainty in the redshift dependence of the es- 
cape fraction of ionising photonsjlnoue et al.|2006||Shapley 
et al | [2006] |Siana et al.1[2010| |Kuhlen fc Faucher-Giguere 
20121. Here, we model the effect of a redshift dependence 

on the evolution of the mass 



of the escape fraction, / csc , 
averaged ionization fraction using three different redshift 
dependences. 



Model 0: /, 
Model I: /< 
Model II: /, 



sc is constant with redshift. 
, sc increases with redshift. 



esc 



decreases with redshift. 



Figure [T] shows the assumed redshift evolution of / csc 
(left panel), together with the resulting mass averaged 
ionization fraction ({xi}) from the different models (right 
panel). Model II completes reionization earlier than the 
other two models, and has the fastest evolution. On the other 
hand, Model I completes reionization later than other two 
models and has the slowest evolution. 

4.2 Halo mass dependence modelling 

It is unknown whether the escape fraction of photons 



Wise & Cen 


2009 


Yajima et al. 


2011 



host halo mass I 


Gnedin et al. 


2008 


|Razoumov & Sommer-Larsen| 


2010 



Here we model the effect of a host 
halo mass dependence of escape fraction on the evolution 
of the mass averaged ionization fraction using two different 
assumptions for the host dark matter halo mass dependence. 

Model III: / csc increases with halo mass. 
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Figure 1. Left: The evolution of / csc as a function of redshift for the three redshift-dependent models 0, I and II. Right: The evolution 
of the mass averaged ionization fraction evolution for the models 0, I and II. 




l°gio("h.,yh-'M ) 



Figure 2. Left: The values of / eS c as a function of host dark matter halo mass for two halo dependent models III and IV. Right: The 
evolution of the mass averaged ionization fraction for the two models III and IV. 



Model IV: / csc decreases with halo mass. 

Figure [2] shows the assumed dependence of / csc on host 
dark matter halo mass (left panel) and the evolution of the 
resulting mass averaged ionization fractions ({xi)). Mod- 
els III and IV show similar evolution of ionization fraction. 
This is because the ionizations are dominated by galaxies 
near lO lo h _1 M , above which galaxies are rare, and below 
which flux is suppressed by feedback processes. The differ- 
ence between the ionization fraction evolution histories of 
the host halo mass dependent models is smaller than for the 
redshift dependent models in Figure. [T] 



4.3 Modelling without supernovae feedback 

For comparison, we also include a model without supernovae 
(SNe) feedback that is otherwise the same as Model 0. This 



model is referred to as NOSN. Note that we adjust / csc for 
the NOSN model as a function of redshift so as to match 
the ionization fraction from Model 0. The resulting values 
of / csc =0.012, 0.017, 0.023, 0.032 and 0.041 correspond to 
redshifts of z ~ 9.278, 8.550, 7.883, 7.272 and 6.712. 



5 THE 21-CM POWER SPECTRUM 

The filtering procedure described in § [3] provides 3- 
dimensional maps of the ionization structure and 21cm in- 
tensity within the Millennium-II Simulation box. From this 
cube we calculate the dimensionless 21-cm power spectrum 

A 2 (fc) = k 3 /(2n 2 )P 21 (k, z)/T (zf (8) 

as a function of spatial frequency k, where P2i(k) is the 
3-dimensional power spectrum of 21cm brightness tempera- 
ture AT(z) (described by eq. (JsJ ) . 
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Figure 3. The evolution of the predicted 21-cm dimensionless power spectra. Colours and line styles show different redshifts as shown 
by the key. 



Figure [3] shows the evolution of the predicted dimen- 
sionless 21-cm power spectra for all models. The predicted 
power spectra show different amplitudes and shapes at the 
same redshift partly because each model has reached a differ- 
ent stage of reionization, and partly because of the relative 
fraction of ionizing photons emitted from large verses small 
galaxies. 

demonstrated that first generation 



Lidz et al. 



([20081 

low frequency arrays like the MWA^] should have sufficient 
sensitivity to measure the amplitude and slope of the 21 cm 
power spectrum. To quantify the effect of assuming an / esc 
that depends on redshift or host dark matter halo mass we 
therefore compare the amplitude and slope of predicted 21- 
cm power spectra for different models. In Figure [1] we plot 
these values as a function of (xi) for central wave numbers 
corresponding to the point on the power spectrum at which 
we evaluate the amplitude and gradient. We choose results 
for k p = 0.2/i _1 Mpc and 0.4/i _1 Mpc, corresponding to the 
range of wave numbers to be probed by the MWA. We plot 
both redshift-dependent (left panel) and host dark matter 
halo mass-dependent (right panel) models for / osc . Interest- 
ingly, models 0-IV which describe variation of / esc with host 
halo mass and redshift do not show significant variation in 
either amplitude or slope. However the NOSN model shows 
significant differences across a range of ionization fraction. 



www.mwa.org 



In particular the amplitude of the NOSN power spectrum is 
lower than for models 0-IV. 



5.1 Observational implications 

Since the ionization fraction is not a direct observable, we 
also plot the progression of a model in the observable plane 
of power spectrum amplitude vs slope. These are shown for 
models 0-IV in Figure [5] again for the two values of cen- 
tral wavenumber k p . The arrows show the direction from 
high to low (xi). To illustrate the potential for detectability 
of this difference we also include error bars for wavenum- 
ber k p =0Ah~ Mpc corresponding to estimates for the MWA 
(specifically an r~ 2 distribution of 500 antennas) (Lidz et al. 



|2008[ ) assuming 1000 hours integration and 6MHz band- 
passes. Figure [5]compares models I-IV and NOSN with the 
Model (which has constant / csc ). As noted in the previ- 
ous section, we find that the effect of SNe feedback strength 
has a much larger influence on the evolution of ionization 
structure during the EoR than the different dependencies of 
escape fractions on redshift or host dark matter halo mass 
considered in this paper. The model with no SNe feedback 
produces curves in the A 2 (fc)-dlnA 2 (fc)/dlnk plane that lie 
at lower values of A 2 (k). Observation of a large A 2 (fe) would 
therefore imply the presence of SNe feedback (or other mech- 
anism which raises the bias of ionizing sources) , irrespective 
of the dependence on ionization fraction within the range of 
models considered. 
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Figure 4. Plots of the evolution in dimensionless 21cm power spectrum amplitude (lower panels) and slope (upper panels) as a function 
of ionization fraction (xi). Predictions are shown for redshift-dependent / csc models (left panel) and halo mass-dependent / csc models 
(right panel). Results are shown for two central wave numbers, k p = 0.4/i _1 Mpc (left) and 0.2h~ 1 Mpc (right) in each case, corresponding 
to the point on the power spectrum at which we evaluate the amplitude and gradient. 



6 CONCLUSION 

Over the next decade we are likely to see the first measure- 
ments of the power spectrum of redshifted 21cm fluctua- 
tions from neutral hydrogen structure during the Epoch of 
Reionization. One goal of these experiments will be to learn 
about the properties of the galaxies that drove the reioniza- 
tion process. It is known that the ionization structure of the 
IGM, and hence the observed 21cm power spectrum, will 
be sensitive to the astrophysical properties of the reionizing 
galaxies. In this paper we have extended the semi-analytic 



model for reionization presented in Kim et al. (20131 to in- 



clude the possibility of redshift or mass dependent escape 
fractions. This model combines the GALFORM galaxy for- 
mation model implemented within the Millennium-II dark 
matter simulation with a semi-numerical scheme to describe 
the resulting ionization structure. We find that an escape 
fraction which varies with galaxy mass and redshift influ- 
ences the structure of reionization. However, we find that the 
effect is smaller than the dominant astrophysical influence 
of SNe feedback. Thus, we conclude that the unknown de- 
pendence on the escape fraction will not influence the ability 
of observations to determine the dominant physics of galaxy 
formation during reionization from the observed 21cm power 
spectrum. 
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Figure 5. Plots of the loci of points in the parameter space of 21-cm power spectrum amplitude and slope. Loci are shown for each 
of redshift-dependent / csc models I and II (top panels), NOSN model (middle panel), and halo mass dependent models III and IV 
(bottom panels) with the Model prediction plotted for comparison in each panel. Results are shown for two central wave numbers, 
k p = 0.2h~ 1 Mpc (top of each panel) and 0Ah~ 1 Mpc (bottom of each panel), corresponding to the point on the power spectrum at which 
we evaluate the amplitude and gradient. The error bars for 0Ah~ 1 Mpc at each point on the Model correspond to estimates for the 
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